Simulating Compact Quantum Electrodynamics with ultracold atoms: 
Probing confinement and nonperturbative effects 
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Recently, there has been much interest in simulating quantum field theory effects of matter and gauge 
fields. In a recent work, a method for simulating compact Quantum Electrodynamics (cQED) using 
Bose-Einstein condensates has been suggested. We suggest an alternative approach, which relies on 
single atoms in an optical lattice, carrying 2/ + 1 internal levels, which converges rapidly to cQED as 
/ increases. That enables the simulation of cQED in 2+1 dimensions in both the weak and the strong 
coupling regimes, hence allowing to probe confinement as well as other nonperturbative effects of 
the theory. We provide an explicit construction for the case / = 1 which is sufficient for simulating 
the effect of confinement between two external static charges. 



Dynamic Gauge theories are in the core of the standard 
model of particle physics, playing the role of the force 
carriers among the matter fields, and therefore are of 
particular significance. It was shown, using lattice gauge 
models and other methods that such gauge theories ex- 
hibit the peculiar phenomenon of confinement of charges 
which is related to non-perturbative effects due to non 
linear interactions in the theory[l -^. Such lattice gauge 
theories are believed to have a non-trivial phase struc- 
ture. The simplest such theory is compact QED (cQED) 
- a U(l) lattice gauge theory, which is believed to mani- 
fest, in 3+1 dimensions, a phase transition between the 
confining phase (for large coupling constant g) and the 
nonconfining Coulomb phase (for small coupling), while 
in 2+1 dimensions it was shown that the theory confines 
also in the weak coupling regime because of nonpertur- 
bative effects [THB]. In non-abelian Yang-Mills theories, 
it is believed that confinement holds for all values of the 
coupling constants. 

Recently there has been much interest in quantum sim- 
ulations of quantum field theories by utilizing methods 
in ultracold atoms and other systems [7 . Models have 
been suggested for simulating dynamical matter fields [8,- 
[TO] , and exotic phenomena manifested by such fields have 
been discussed [TTI-FI3] . However, less progress has been 
achieved for dynamic abelian and non-abelian gauge the- 
ories. Dynamic gauge theories involving spin-half states 
have been discussed in [14|(15J. Since the electric fields 
in such models can obtain only two values, such mod- 
els are unable to manifest the effect of electric flux-tubes 
(but rather of different "strings"). Coloumb phase sim- 
ulations have been suggested with molecular states [16^ , 
and BECs in optical lattices [17 . 
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In a recent work [TS^, we have obtained, using Bose 
Einstein condensates in an optical lattice, an effective 
theory of a dynamic U(l) gauge theory, manifesting con- 
finement of external static charges, with observable elec- 
tric flux-tubes. In the present paper we suggest an alter- 
native approach for simulating gauge theories in terms of 
a Spin-Gauge Hamiltonian Hi (defined in equation [T]), 
which describes interacting single atoms with internal 
levels playing the role of angular momentum multiplet 
— / < m < instead of BEGs as [18 . We will show that 
for large values of / the Spin-Gauge Hamiltonian rapidly 
converges to the standard abelian Kogut-Susskind model 
of cQED for both the weak and strong coupling 
regimes. Hence this model is able to simulate the effect 
of confinement as well as non-perturbative effects in the 
weak coupling regime, which give rise to it in 2+ ID. As 
a first step in realizing the models Hi we shall propose a 
method for constructing the case / = 1 which is sufficient 
for observing the effect of confinement of static charges 
in the strong regime. First, using the methods of [20 , we 
construct an effective "generalized XXZ" model. Then, 
adapting the ideas of [18 , the gauge invariance is intro- 
duced to the system with a second effective calculation. 
The single atoms implementation might be easier exper- 
imentally than the BEG approach, as it does not rely 
on the overlaps of local single-particle wavefunctions and 
thus the energy scales in the Hamiltonian may be larger. 
Moreover, only three atomic species are required, and 
they can populate every link, unlike in the BEG model. 

Let us consider a 2D square lattice with single atoms 
which carry 21+1 internal states, located on the links and 
described by the Spin-Gauge (SG) Hamiltonian: 
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where n are the lattice's vertices, /c = 1,2 are the lat- 
tice directions, whose corresponding unit vectors are 1,2. 
For example, L^^ is the z component of the spin on 
the link emanating from the vertex n in the kth direc- 
tion (The generalization to a 3D lattice is straightfor- 
ward). ^ is a constant. This should be compared to 
the abelian Kogut-Susskind Hamiltonian |2l [19] Hks 
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like the Kogut-Susskind Hamiltonian, in our case we are 
dealing with 3D angular momentum operators. Never- 
theless, for large values of the first quadratic term in 
the Hamiltonian conincides with the electric part of the 
abelian Kogut-Susskind Hamiltonian , with the z com- 
ponents of the angular momentum playing the role of an 
electric field and the second, quartic part with the 
magnetic part of the Kogut-Susskind Hamiltonian. This 
can be seen qualitatively when considering the matrix el- 
ements of \m\ <C for which 1^, ^) ^ |^ ^ i l)^ 

similarly to e^*"^ |m) = |m ± 1) in the Kogut-Susskind 
model. We shall test this equivalence quantitatively for 
the case of a single plaquette (figure [T]) . 

For every /, the Spin- Gauge Hamiltonian manifests 
a local U(l) gauge symmetry, that is generated by the 
local operators (defined on the vertices of the lattice) 
Gn = E^/c^L (where A^/n /^^k " /") which 

k 

commute with the Hamiltonian: for a given vertex, Gn 
trivially commutes with all the plaquettes which do not 
contain n. As for the other four plaquettes, the com- 
mutation relation is zero, since [Lz^L±] = 

For 
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example, Gn, ^V,n^+,n+i^i,n+2^^ 
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Static external 



charges |{Qn}) are introduced to 
the system by fixing a subspace by the constraint 

Gn|{Qn})=gn|{gn}). 

In order to have something useful for simulations, we 
would like to have that Hi Hks for large /s, suffi- 
ciently fast. Thus we shall consider a comparison be- 
tween the Spin-Gauge Hamiltonian with a constant / and 
a truncated version of the Kogut-Susskind Hamiltonian 
with —1<E<1. It is straightforward to see, using 
perturbation theory in ^ that in the strong limit of 
the Hamiltonian {g ^ 1) the ground states of the Spin- 
Gauge and Kogut-Susskind Hamiltonians coincide up to 
a certain order in the perturbative expansion, depending 
on / and the charge distribution. On the other hand, 
in the weak coupling limit, we shall examine the effect 
of truncation in a nonperturbative manner, for a single 
plaquette system. 

Case of a single plaquette. Consider a single plaque- 
tte with two opposite unit static charges in the lower 
vertices (see figure [l]^). Using the gauge invariance and 
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FIG. 1. Single Plaquette Plots, (a) Graphs of {Ei (g)). Black 
- the calculation of [5] - small coupling approximation for reg- 
ular abelian Kogut-Susskind Theory. Blue - exact calculation 
for the truncated theory, for for / = 1,...,20. Red - Exact 
calculation for the Spin-Gauge theory, for / = 1, 20. It can 
be seen that the curves start to merge for a small g and /. 
The value 1 refers to the flux-tube and | to the longitudal 
part, (b) Graphs of Pi (g), for / = 1,...,20. It can be seen 
that even for small values of /, Pi (g) approaches 1 for a small 
g. (c) The one plaquette system we use in the demonstration. 



Gauss's law, a possible gauge- invariant basis of states 
is |m) = |m, m — 1, 1 — m, 1 — m) (m = — / + 1,...,/) 
where these are the eigenvalues of the electric field on 
each link, from the lower one, counter-clockwise. Re- 
lying upon the results of Drell et. al. in [5], the 
ground state of this system, for weak coupling, is given 
by a Bloch-like wavefunction in the tight-binding limit 

X{0)= E 



27rn) ^ Q magnetic 



field in the plaquette. On the link connecting the charges, 

(El) = 1^^ (^^) ™^ understood as 

the contributions of two parts: One is the trivial contri- 
bution of the longitudal, classical, static Coulomb field 
(|). The second part is much more interesting: it is non- 
analytic in ^ = and therefore it is nonperturbative in g. 
It is this type of mechanism which is responsible for con- 
finement in the weak regime in large 2D spatial lattices. 
In figure [l^, (Ei) is plotted as a function of ^ - both 
for the truncated Kogut-Susskind Hamiltonian and the 
Spin-Gauge Hamiltonian for several values of and the 
analytic result of [5 as well. For ^ < 1, the truncated 
and Spin-Gauge results conincide with each other and 
with the analytical approximation for / ^ 2,3. For large 
values of ^, the ffux-tube value of 1 is reached already for 
/ = 1. 

In order to understand the effect of truncation, we ex- 

oo 

pand this state in the m basis - |x) = |^) = 

m= — oo 

oo 2 / 3 \2 

^ ^-9 [^-4) \rn), and take only + 1 < m < /. 
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FIG. 2. (a) The simulation lattice. As explained in the text, 
the atoms (yellow circles) are aligned along the links of a 
2d square lattice, with basis vectors 1 and 2. On the blue 
lines the tunneling rates between neighboring atoms are (in 
Hamiltonian (2)), and on the red lines - t^. In Hamiltonian 
(3), z — z interaction are between nearest neighbors along 
both red and blue lines, and x — x, y — y interactions are 
only along the red ones, (b) The flux tube (zeroth order 
in perturbative expansion, in the strong limit of the Kogut- 
Susskind Hamiltonian) connecting two opposite charges, (c) 
The flux tube in the language of our simulating system. 



late the probablity to be in the truncated subspace: 



Pi{9) = 



(xlx) 



= -^ + 1 



as one can see in figure 



l3, 



this function approaches 1 quickly even for small finite 
/s, which means that the truncated theory still shows the 
same effect for small ^s. 

Simulating the I = 1 Spin-Gauge Hamiltonian. Let us 
consider a 2d square optical lattice[21 , whose minima 
conincide with the links of the square lattice of the Spin- 
Gauge Hamiltonian (the generalization to 3d is straight- 
forward). Each minimum is populated by a single atom 
of three different atomic levels, forming a / = 1 spinor 
(see figure [2|l). We use it to develop an effective theory 
[22] . which manifests confining flux tubes similar to the 
ones in cQED. We first turn to the derivation of an effec- 
tive generalized spin-1 XXZ Hamiltonian on this lattice. 

The Hamiltonian describing the dynamics of three 
atomic species a G {+,0,—} on such an optical lattice 
contains same-species tunneling terms along the 1,2 di- 
rections and the diagonal directions as well, and on-site 
spin and number dependent terms. A pictorial represen- 
tation of the interactions can be found in figure [2^, and 
a more deatiled description of the lattice structure and 
the optical potential can be found in the supplementary 
material. Generalizing from [20 we obtain the Hamilto- 



As a measure of the accuracy of truncation we calcu- nian: 



+ E + yE^n {N^ - 1) + yE {^"^ - 2^n) 

n,k,a n,k n,k 



where for horizontal links the annihilation operators are 
an,a, and the number operators - a — a^n,^, and 
for vertical links -b^^a and ^ = bl^^bn^a] = 
^^na total ou-sitc spiu (scc [20]). 

a 

We set the parameters ^ = |+2A + /i^ 
2A ((7n + ^n+k) , Ao = 0, with < /7o, /72 , /i < A < 
t/o, ^2, and are integer C-numbers which will be later 
related to the static charges. We also introduce a new 
variable satisfying U2 = zUq. 

Derivation of a generalized XXZ Hamiltonian. The 
first effective calculation, which leads to a generalized 
XXZ model, is similar to the one in [20 . Since the U 
local terms are much larger than the others, it is reason- 
able to obtain, perturbatively, an effective Hamiltonian 
around them [22 . Unlike in [20 , we do not include the 
small local terms ^ a^n. a constraining part 

n,/c,Q! 

of the Hamiltonian, but rather treat them as the first or- 



der contribution to the effective Hamiltonian. We divide 
that into two parts: one is the A,/i,(7n dependent part, 
which we put aside at the moment; The other 5 depen- 
dent part will be used to construct "two-site" connected 
local diagonal Hamiltonians: the energy contribution of 
this part from each link will be equally distributed among 
the six connections it has with other links (to whom it is 
connected by tunneling). 

The tunneling rates are chosen to be real: = 

fs _ +d _ +d _ 1 / £/o(24A-5e)(e+24A) +s _ n +d _ 
t_ — t+ — C_ — 4 Y , Iq — — 

^/iiS^' where g « e « 1^ « A and 
t^ < 4- Finally we set S = -12A - e and z = | - ^. 
Then we can apply the effective calculation as in [20] . 
Collecting these resutls with the A, /i, dependent part, 
we get, up to a constant energy, the first effective Hamil- 
tonian, 
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^e// = 2^ E L,,nL,,r.'+nJ2{L.,nL.,n'+Ly,r.Ly,r.') + ^[{2X + ll){Ll^f -2\{^^ (3) 

str-\-diag diag n,k 



where the z — z interactions are between hnks who share 
a vertex and the x — x^y — y only between hnks who share 
both a vertex and a plaquette (see figure [2|i). This is an 
example of a frustrated XXZ model Hamiltonian, which 
is of interest of its own. One can easily check that the 
scale hierarchy is not violated. 



which is similar to the Hamiltonian obtained by us pre- 
viously [18] [23] and from which, due to the scale hi- 
erarchy A ^ ^, we shall obtain an effective Hamil- 
tonian as in [18 . The constraint will be Hq (Gauss's 
law), and its ground sector contains the states of rele- 
vance for us. He commutes with it and hence becomes 
the first order of the effective Hamiltonian, and from 
Hr we get two conributions: one is the gauge invariant 

plaquette term Hb = -^Y. , 2^"^n^n+i + H.C.) . 

m,n ^ ^ 

The other one is due to the finiteness of the angular mo- 
mentum representation matrices, but it is diagonal (and 
hence gauge invariant) and thus introduces a negligible 
first order correction to the energy but does not change 
the zeroth order ground state (the fiux tube): H'^ = 
E (1+) (+1 + |0) {0|)„ ® (|0) (0| + I-) (-!)„,. Note 

diag 

that as / ^ oo (in a truncated Kogut-Susskind theory), 
these terms approach identity matrices and act as an ig- 
norable costant energy, and hence this term did not ap- 
pear in the infinite case. An example for the emergence 
of gauge invariance as the constraint gets stronger is il- 
lustrated in figure |3] 

Next we perform the change of variables [18]: 

^^,n ^ i-'^T'^^'^Lin^ <?^n ^ (-1)^^+^^^^ and Qn = 

(^_2)ni+n2^^ (which also swaps the L± operators on odd 
vertices). This change of signs is needed in order to get 
the correct signs for the Gauss's Law constraint and the 
interactions of the Spin-Gauge Hamiltonian. Plugging it 
into the effective Hamiltonian we get 

^e//=A^E(^U'+^B + ^B (5) 

and if we rescale the energy, using a = ^-/i = "^^^ ^ , 

we get the that a~^H^^j is the / = 1 Spin-Gauge Hamil- 
tonian equation ([l]) plus the irrelevant term of a~^H'^. 
Hence, because of the sign inversions between the Spin- 
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Imposing gauge invariance on the system. In the sec- 
ond step, we employ the method of ^8j. Define Cn = 

LF' and rewrite the effective Hamiltonian as 



(4) 
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FIG. 3. An example of the emergence of gauge invariance in 
the second effective calculations for a single plaquette system 
(/=1). Two unit charges are located in the lower vertices. The 
parameters are chosen such that for a large A the system will 
be in the strong limit: fi — 1, Q — 0.1. The expectation values 
of Lz along two links, as well as of the gauge transformation 
generators G at two vertices are plotted as a function of A for 
the ground state. It can be seen that as A grows, G converge to 
the local charges and hence Gauss's law and gauge invariance 
are introduced, and Lz approach the values of the lowest flux- 
tube state in the strong limit. Note the discussion on the signs 
in the text. 



Gauge and the simulating models, the manifestation of 
an electric flux-tube between two conflned static charges, 
to zeroth order, will be a line of alternating +, — states 
of the atoms along the links connecting the two charges 
(see flgure[2]3). Some speciflc suggestions for initial state 
preparation and possible measurements are presented in 
the supplementary material. 

In this paper, we have presented a new method to sim- 
ulate cQED using ultracold atoms in optical lattices. We 
believe that this method may be experimentally acces- 
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sible in the near future. Although we have constructed 
a reahzation of Hi, which ahows for simulation of con- 
finement around the strong coupling limit, the rapid con- 
vergence of Hi to the Kogut-Susskind model suggests an 
avenue towards the simulation of the nonperturbative ef- 
fects of the weak coupling limit as well as phase transi- 
tions in 3+1 dimensions. It would be intriguing to study 
the inclusion of dynamic matter fields in the model, which 
would lead to a full simulation of cQED. 

Acknowledgments. The authors would like to thank L. 
Mazza for helpful discussions. BR acknowledges the sup- 
port of the Israel Science Foundation, the German-Israeli 
Foundation, and the European Commission (PICC). IC 
is partially supported by the EU project AQUTE. 

SUPPLEMENTARY MATERIAL 
The superlattice structure 

In the paper, we have suggested a method for simu- 
lation of the / = 1 spin-gauge hamiltonian, as a simu- 
lation model for cQED. The effective derivation of the 
gauge- invariant Hamiltonian relied on Hamiltonian (2) 
in the paper, which is a generalization of the model of 
[20 , where in section IV and appendix A the method of 
obtaining the tunneling rates is explained in detail. 

In our case, the system is 2d and thus the generaliza- 
tion should be explained. As in [20 , the different hop- 
ping rates are achieved using a suplerlattice structure, 
in which unoccupied auxiliary minima, which are cou- 
pled by lasers to the "regular" occupied atomic minima, 
are adiabatically eliminated forming the desired hopping 
processes. 




FIG. 4. The Superlattice Structure A rotated plot of the op- 
tical lattice (45° degrees rotation from figure 1 in the paper). 
Yellow circles denote the occupied sites (which are on the links 
of the paper's model, which are the blue lines, and on the ver- 
tices here). Blue and Red circles denote the auxiliary sites of 
two different depths, through which tunneling processes be- 
tween yellow sites are generated. 



In our model, the atomic minima of interest are located 
on the links of a spatial 2d lattice (see figure 2a in the 
paper). In figure 1 of the supplementary material, we 
present a rotated version of the lattice, which includes 
the auxiliary sites as well. As one can see in the figure, 
there are two types of auxiliary sites, with two different 
depths. This can be achieved using an optical potential 
of the form 

V (x) = -Vo[cos^ (2kx) cos^ (2ky) + 
cos^ (kx) cos^ (2A:?/) + 
cos^ {2kx) cos^ {ky) + 




where Vb, > (see figure 2). 



FIG. 5. A Contour Plot of the optical potential. The minima 
are colored according to figure 1. The yellow minima are the 
deepest. 

To achieve the requested experimental scheme, one 
laser connects (along the blue lines in figure 1) the 
m = ±1 levels of "real", "yellow" sites and the same levels 
of the neighboring "blue" auxiliary sites, giving rise (us- 
ing adiabatic elimination) to the (equal) tunneling rates 
tj^^t^ between "yellow" atoms. Another laser connects 
(along the red lines in figure 1) the m = levels of "yel- 
low" sites with the same levels of the "red" sites. These 
lasers are responsible to t^ in a similar manner. No other 
lasers are required, since we need = 0. 

State preparation and possible measurements 

Suppose, first, that the bosonic interactions {Hb^H'^) 
are turned off (17 = O = 0) - i.e., the system is in 
the strong coupling limit. The gauge field vacuum (when 
for every vertex n, = 0, does not contain any excited 
links: all the atoms are in their m = state. This is the 
exact QED vacuum in the extreme strong limit, since no 
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plaquette terms are present. Thus the system has to be 
initiahy prepared in this state, which is an eigenstate of 
Hq + He (and hence they do not need to be turned on 
adiabaticahy - if they were initiahy of, any way of turning 
them on can not change the state, but the Gauss's law 
constraint wih be imposed). 

After obtaining this state as a first stage, one can either 
do nothing or obtain some interesting initial states in the 
following methods; Most of them involve an adiabatic 
increase of 1^, which is an adiabatic increase of t^'. 

1. Remain in the strong limit; Do nothing else. 

2. Obtain the QED vacuum (up to first order): in- 
crease adiabaticahy, keeping ^ <C /i. This 
would effectively mean an adiabatic increase of Hb 
(and i^^), dressing the system to the desired state. 
Then, if Vt remains fixed after the adiabatic evo- 
lution, QED dynamics (to first order) will govern. 
One can use different Vts throughout the lattice, 
thus "focusing" the loops to certain areas (for ex- 
ample, to the center). In this case, the corrections 
to the strong limit vacuum can be calculated per- 
turbatively. It can be shown using the theorem in 
the next section, that for / = 1, the perturbative 
series of the vacuum states for the Kogut-Susskind 
and the Spin- Gauge theories are identical to first 
order; The first order correction is just a superposi- 
tion of single excited plaquettes all over the lattice. 

3. Obtain a "loop sea" by increasing Vt adiabaticahy 
again, but to ^ ^ ji. Thus the system will be 
dressed to a state with lots loops, but with larger 
amplitudes. Then, when the state is ready, lower 
the value of Vt to obtain strong regime QED dy- 
namics. 

4. Create "Manually" zeroth order excitations: us- 
ing single-addresing lasers |25], one can cre- 
ate closed flux-loops, or, if there are nonzero g'nS, 
one can generate the required zeroth order flux- 
tubes, which are made of alternating signs atoms, 
as explained in the paper. Afterwards Hb can be 
turned on adiabaticahy to obtain QED dynamics 
(note that in this case the highest accurate order is 
the zeroth order, as explained in the next section, 
but some corrections to it do appear). 

A qualitative glimpse on the structure of the state can 
be achieved if the density of atoms of each species is 
measured and " frozen" (by switching off ~ after the 
measurement). The state collapses to a single element of 
the superposition, and fiux tubes and loops be observed. 
Being in the confining phase, only finite flux lines can be 
seen: either closed loops (single plaquettes, since this is 
the leading order contribution to the state) or flux tubes 
between static charges, if there are any. In the latter case 
the flux-tube prepared in zeroth order should remain, but 



with an additional loop somewhere around the lattice. 
Although the probability for a particular loop is small, 
as a perturbative correction, if the system is large enough, 
the probability for a loop (somewhere around the lattice) 
becomes larger. This is unless the initial state was the 
"loop sea" - in that case, other interesting loop structures 
may be observed. 

One should also note that the predictions given are 
based on perturbation theory considerations, and thus a 
realization of these measurements may serve as an inter- 
esting experimental check of them. 



Accuracy of Hi eigenstates in perturbation theory 

Theorem. Let {N^} be operators with eigenvalues 
n/e G Z and eigenstates \nk) satisfying Nk\nk) = nk\nk) • 
Let |a/e, be lowering and raising operators, such that 

ak\nk) oc \nk - I) ,al\nk) oc |n/c + 1). Consider the 
Hamiltonian 

H = Ho{{Nk}) + Xv(^{ak,al}) (7) 

where A <C 1, and the maximal power of any lower- 
ing\raising operator in F is m, and a truncated version 
of the Hamiltonian, 



H = Ho [{Nk]) + W {^[ak,ai}) 



(8) 



where the eigenvalues and eigenstates of |^/e| satisfy 

—l<nk<l and ak — 0,a^ |/) = |0). Suppose that 
in the ground state [26 of i^o ({A^fc}), the highest 

eigenvalue of an {A/'^} operator (in an absolute value) is 
no. 

G ) of H., H respec- 



(9) 



Expand the ground states \G) 
tively in time independent perturbation theory. 



oo 

\G) = |g(°)^ + ^a"|g(")^ 

n=l 



G = 



n=l 



then for every n < 



(10) 



(11) 



Proof. Start with n = 0. In order to have |G^^^) = 
given that |G^^^) contains the maximal eigen- 
value no, we need that I > uq. Next, for n = 1, the high- 
est possible maximal eigenvalue (abs. val.) is no + m, 
by applying the raising\lowering operators in a maximal 
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power of m. Thus one needs, for first order accuracy, 
that / > no + m. Continuing to higher orders, we get 
that for nth order accuracy, one must have / > no + nm, 
and hence, finally, n < QED. 

For the Hamiltonian of compact lattice QED, the theo- 
rem applies to the strong coupling limit. There, Hq is the 
electric part, and V is the magnetic part. Since the pla- 
quettes contain a single power of raising\lowering opera- 
tors, m = 1. In the absence of external charges, no = 0, 
and thus n < i.e., if the electric fields are truncated in 
a value of this will be the highest accurate order. If 
there are, for example, two opposite static charges con- 
nected by a fiux tube, no = 1, and thus n < / — 1, i.e., to 
be accurate to order n, the truncation limit must be at 
least n + 1. 

For example, approximating the system with spin-1 
spinors, one gets accuracy to first order in the pure gauge 
case but only to the zeroth trivial order in the case of 
external charges, for which an approximation using / = 2 
is needed in order to get first order accuracy. 

For the sake of a complete discussion, one should note 
that in order to apply the theorem to our class of trun- 
cated operators, we should also consider the angular mo- 
mentum nature of the spin-gauge model. In the / = 1 
case, although the ladder operators are not unitary, their 
nonvanishing matrix elements are identical and hence 
this problem can be disregarded. We shall coment that a 
"cure" for this problem (using some local counter terms) 
exists for / = 2. 
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